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Abstract. The local moment approach is extended to the orbitally-degenerate 
(SU {2N)) Anderson impurity model (AIM). Single-particle dynamics are obtained 
over the full range of energy scales, focussing here on particle-hole symmetry in the 
strongly correlated regime where the onsite Coulomb interaction leads to many- 
body Kondo physics with entangled spin and orbital degrees of freedom. The 
approach captures many-body broadening of the Hubbard satellites, recovers the 
correct exponential vanishing of the Kondo scale for all A*^, and its universal scaling 
spectra are found to be in very good agreement with numerical renormalization 
group (NRG) results. In particular the high-frequency logarithmic decays of 
the scaling spectra, obtained here in closed form for arbitrary A'^, coincide 
essentially perfectly with available numerics from the NRG. A particular case 
of an anisotropic Coulomb interaction, in which the model represents a system 
of A'^ 'capacitively-coupled' SU (2) AIMs, is also discussed. Here the model is 
generally characterised by two low-energy scales, the crossover between which is 
seen directly in its dynamics. 
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1. Introduction 

The Anderson impurity model (AIM) [1] plays a central role in our understanding 
of correlated-electron physics. It has long been used to explain the properties of a 
diverse range of systems where strong, local Coulomb interactions arc paramount [1], 
most recently in the area of electronic transport through nanostructures such as 
semiconducting quantum dots, carbon nanotubes and single-molecule devices [2,3]. 

The original AIM [4] describes a system of spin-i electrons, tunnelling between 
a metallic conduction band and a single impurity level on which they experience a 
local Coulomb repulsion This interaction renders the problem highly non-trivial, 
but the physics of the model has been uncovered over several decades by a range of 
techniques, see e.g. [1,5-19]. Most significantly, the Kondo effect — the quenching of 
the impurity spin via formation of a complex many-body state between the impurity 
and conduction band - is now rather well understood. 

In this paper, we consider the SU {2N) generalisation of the AIM [1] , in which 
the spin-^ electrons of the original model are replaced by particles with 2 A'' internal 
'flavour' degrees of freedom. As well as being important from a theoretical point of 
view, the N = 2 case in particular has direct relevance to recent experiments on carbon 
nanotube quantum dots [20-24]. The flavour index of the particles here corresponds 
to a combination of electron spin and orbital indices, and the model exhibits an S't/(4) 
Kondo effect in which spin and orbital degrees of freedom are entangled [25,26]. 

The SU (2A/') AIM has appeared in the literature in many guises over the years. 
In the f/ — > oo limit, where the impurity can be at most singly-occupied, the model 
becomes somewhat simpler to analyze. It is integrablc, and hence can be studied 
using the Bethe ansatz [9,27-29]. Various approximations falling under the umbrella 
of 'large-iV expansions' [10-16] have also been developed. These become exact in the 
limit N — > oo, where their thermodynamics agree with the Bethe ansatz results [1]. 
Generalising large- techniques to finite-[/ is however complicated [30-35], rendering 
this still an active area of research today [36] . Recent progress has also been made in a 
different direction using a slave-rotor mean-field theory [37] , which recovers the correct 
exponential form of the Kondo scale for arbitrary N at the particle-hole symmetric 
limit. 

In addition to analytical approaches, SU{2N) models have also been studied 
numerically. Thermodynamics and dynamics have been calculated using quantum 
Monte Carlo (QMC), sec e.g. [38]. While QMC can in principle handle any degeneracy 
A^, interaction strength U and temperature T [38], the method is computationally 
rather expensive, and cannot reach the exponentially-small low energy scales on 
which the Kondo effect is fully manifest. The numerical renormalization group 
(NRG) [7,8,39] has also of course been applied with great success to the SU{2) model 
(for detailed reviews see e.g. [1,39]), as well as to finite-C/ SU{4) AIMs, e.g. [24,40-43]. 
NRG is numerically exact down to the lowest energy scales, although the scaling with 
N of its CPU and computer memory usage is such that at present it is realistically 
limited in practice to the SU{2) and SU{A) cases. 

Here we analyze the SU(2N) model using the 'local moment approach' (LMA), 
considering explicitly the particle-hole symmetric limit where the impurity is occupied 
by N particles. The LMA - developed originally for the SU{2) AIM [19,44-49] and 
subsequently extended to encompass both the pseudogap [50-53] and gapped [54, 55] 
AIMs, as well as lattice-based models within dynamical mean-field theory [56-61] - 
has been shown to circumvent many of the traditional deficiencies of previous many- 
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body approaches. It can handle all interaction strengths, from the non-interacting 

limit up to the 'strong-coupling' regime where electron correlations dominate; and 
recovers the correct physics at both high and low energies, the latter being the Kondo 
effect described above and its ultimate Fermi-liquid description on the lowest energy 
scales. It is also physically transparent and computationally straightforward; many of 
its results can in fact be obtained analytically, including non-trivial closed formulae 
for single-particle dynamics of the model in the strong-coupling regime. 

We explain in the following that generalisation of the LMA to the SU{2N) AIM is 
quite natural (and does not require significant additional computation when applied 
to arbitrary N). In the strong-coupling regime, N essentially appears only as an 
additional prefactor to the LMA self-energies, and hence many of the known LMA 
results for the SU{2) model can be extended quite straightforwardly. In particular, 
we show that the LMA recovers the correct exponential dependence of the Kondo scale 
on N, and derive a closed expression for the 'tails' of the imiversal SU{2N) single- 
particle scaling spectrum that agrees excellently with NRG results for the SU{2) and 
SU{4:) cases. 

The paper is structured as follows. Section 2 defines the SU{2N) AIM, shows 
its connection to a model of N capacitively-coupled SU{2) AIMs, and defines the 
particle-hole symmetric case on which we focus here. In section 3 we introduce the 
LMA, starting with an examination of the non-interacting and 'atomic' limits, before 
moving on to the mean-field level of unrestricted Hartree-Fock and then the LMA 
itself. Section 4 contains a summary of our results for the model, both numerical 
and analytic as mentioned above, together with direct comparison to NRG results 
for single-particle dynamics in the SU{4) and SU(2) cases. We end the paper by 
discussing briefly the physics of the model away from SU{2N) symmetry [section 5], 
and provide a short conclusion in section 6. 

2. Model 

In conventional notation the Hamiltonian for the SU{2N) Anderson impurity model 
is [1] 

where = fj^f^, the flavour index m takes on 2A'' distinct values, and the prime 
in the third sum means m' = m is excluded. The Hamiltonian thus describes a non- 
interacting host band of 'flavourful' particles (first term), plus a localised impurity 
orbital with onsite Coulomb repulsion U (second and third terms), with tunnel 
couplings between the two (final term). 

Equation (1) highlights clearly the inherent SU{2N) symmetry of the model. 
But it will also be instructive to rewrite the Hamiltonian in a way that reveals its 
connection to a problem of spinful SU (2) electrons. To this end, we regard m as a 
composite index, m = {i,a), where i takes on N integral 'site' indices from 1 to 
and cr = t or I represents an SU{2) spin index. Sums over m then become joint sums 
over i and a, and (1) can be rewritten as 

H=Y1 'Aia^ua + e E + 1 E' + ^ E {"iiJicr + ii-c-) (2) 

k,i,(T i,(T i,j,a,a' k,i,(T 
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and a = a' is excluded. The interaction term can be simplified by writing 



with rijo- = fj^f^^, and the prime on the sum means that the term with both i = j 



E -EE^'. + EE^^v (3) 

i,j,a,a' i,j a, a' i^j a,a' 



(with Sab the Kronecker delta), such that 



H 



E' 

k,(T k,CT 



where fii = J2a- ^^^^ form, the Hamiltonian describes A'' equivalent capacitively- 

coupled SU(2) Anderson models, in the limit of identical onsite and intersite Coulomb 
interactions. The LMA developed here can in fact be straightforwardly generalised to 
the case when the intersite interaction is less than U ; we comment briefly on the main 
results of such in section 5, but until then focus on the fully-symmetric SU{2N) limit 
corresponding to (1). 

For relative simplicity we also consider the model in the particle-hole symmetric 
limit, with a standard flat, infinitely- wide conduction band with density of states po [!]■ 
The latter automatically ensures particle-hole symmetry of the host; the impurity part 
of (1) under a particle- hole transformation — > — becomes 



^[-e-(2iV-l)C/]n„-h-^'n„n™, (5) 

and thus the model is fully particle-hole symmetric when 

e = -(iV-i)C/. (6) 

All es that appear in the following are implicitly tied to U in this way. For later use, 
note that the Hamiltonian (1) in this limit reduces (modulo an irrelevant constant) to 
the capacitive charging form 

^ = E 'Am'h.m + - A^)^ + F E (cL/„ + h-C.) , (7) 

k,m k,m 

with N = J2m total impurity charge operator. We add here that the 

subsequent approach is not in fact restricted to the particle-hole symmetric limit 
per se; the LMA has already been used successfully to describe the SU{2) AIM away 
from particle- hole symmetry [44], and generalisation of the approach below to the 
particle-hole asymmetric SU{2N) model can be developed in a similar vein. 

To examine the dynamics of the SU{2N) AIM, the central quantity of interest 
is the time-ordered, flavour-m impurity Green function Gm{i^), defined as usual 
[62] by the Fourier transform of G„(f) = -i(T{./„(i)//„ (0)}) (with f the Wick 
time-ordering operator), and with corresponding single-particle spectrum Dm{<^) = 
—w~^ sgn(a;) ImGm(t^)- It is this on which we naturally focus in the paper. 



3. Local moment approach 

Before we develop the LMA for the SU{2N) AIM, it is instructive to comment briefly 
on two special cases: the non-interacting limit of U = (but arbitrary V), and the 
atomic limit with V = (but arbitrary U) . While only glimpses of the physics of the 
full model are seen in these rather trivial limits, our main piirpose in considering them 
is that any credible approximate theory for the model should recover both limits as 
particular cases. 
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3.1. Non-interacting limit 

In the non-interacting hmit, U = 0, the impurity Green function (here denoted by 
9m{'^)) is readily obtained in closed form, in direct analogy to the SU{2) case [1], viz. 

9m{i^) = . \. ^ (8) 

with Lv~^ = ui + 10+ sgn(ti;). Here A(a;) = — Ck) is the host-impurity 

hybridisation function, reducing to A(w) = — isgn(w)Ao for the conventional wide- 
band host considered explicitly, where Aq = Trpol^P (with po the conduction band 
density of states, as above). The spectrum of (8) is simply 

A^ 

7rAorf^(a;)= ° ; (9) 

the Lorentzian form of which reflects the lifetime broadening (of width Ag) of the 
impurity level due to its hybridization with the host band. This introduces a natural 
energy scale to the problem, and it is convenient to use Aq as the unit of energy in 
the rest of the paper (unless stated otherwise). To this end, we define the reduced 
Coulomb interaction U = U/Aq and the reduced frequency u) = cu/Aq. 

3.2. Atomic limit 

The atomic limit {V = 0) is trivially soluble for any value of U, since the impurity 

number operator A'^ commutes with H. The basic physics follows directly from the 
Hamiltonian in (7). When U > 0, the ground state is necessarily degenerate. The 
impurity is occupied by N particles in its ground state, and these can be distributed 
amongst the 2N orbitals in Wn = (^^) distinct ways. Single-particle excitations 
connect the degenerate ground state manifold to further N ± 1-particle degenerate 
manifolds, all lying {7/2 higher in energy. 

In the absence of any special preparation, the particular ground state adopted 
by the system is equally likely to be any one of the Wn possibilities. The single- 
particle excitation spectrum for a given flavour m is then the statistical average of 
the individual flavour-m spectra of all the Wn ground states. In precisely Wn/2 of 
these, there will be a particle in flavour orbital m and hence:}: a 'removal' pole in the 
spectrum at u) = —U/2. In the other Wn/2 cases the mth flavour orbital will be 
unoccupied, leading to an 'addition' pole at w = +U/2. As a result, the spectrum in 
the atomic limit consists of two poles of weight ^ lying at a; = ±U /2, i.e. 



„ / N v=o 1 

Gm{uJ) = - 



1 1 

+ 



(10) 



independently of N. 



3.3. Mean field 



Next we turn to an analysis at static mean-field (MF) level. By analogy to unrestricted 
Hartree-Fock (UHF) for the SU{2) case [4], we allow for the most general form of local 
symmetry-breaking, writing rim = {fim) + for all flavour indices m. Substituting 

I All energies are measured with respect to the Fermi level. 
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this into (1) and neglecting the fluctuation terms quadratic in Snm^fim' gives, to within 
an irrelevant constant, 

^ = Y1 ^kcLck„ + J2 ^mflm + ^ X] + ^-c) , (11) 

with 

em = e + U X {fim-). (12) 

Under a particIe-hoIe transformation N 2N — N, and since H is invariant under 
this transformation (see (7)), particle-hole symmetry implies the single constraint 

Y.{nm)=N, (13) 

m 

SO that there arc 2N — 1 possible channels in which symmetry can be broken at MF 
level. Use of (13) in (12), along with (6), yields an expression for em in terms of (fim) 
alone, viz. 

em = U{^-{hm)). (14) 

Since the Hamiltonian in (11) is separable for each flavour m, the MF impurity 
Green functions follow in direct analogy to (8) as 

= -T — ■ (15) 

L0+ - em + i sgn(a;) Ao 
The mean number of flavour-m particles on the impurity is then given by integrating 
the spectrum Vmit^) = —Tr~^ sgn(a;) Im^m(i^) up to the Fermi level: 

(16) 



{fim) = J 'Dmiuj) dw = ^ - tan ^ (^^^ • 



Together with (14), this implies the following MF self-consistency equation for each 

subject to the constraint of (13). 

The simplest solution of (17) and (13) arises when {hm) = \ for all to. This is 
restricted Hartree-Fock (RHF), in which the SU{2N) symmetry of the impurity is left 
entirely unbroken. In this case, by (12), all the CmS vanish, and hence (15) reduces to 
the non-interacting Green function gm{u}) [(8)] for all to. But while RHF is a (trivial) 
self-consistent solution of the MF equations, it is not the solution of lowest energy. 
By breaking the symmetry of (1), variationally lower solutions can be generated [4]. 

We have analysed the broken-symmetry MF solutions in detail. As now shown, 
the self-consistent solutions have the 2N mean-field energies partitioned into two 
groups of A'' degenerate levels. We label the level energies of the two groups as e+ 
and e_, and the mean occupation of each level as {fi+) and (n-), respectively. For 
convenience in what follows, we simply refer to the two groups of N levels as the 
and '— ' groups, respectively. 

Upon partitioning the levels between these two groups, the 2A'' equations in (17) 
reduce to just two: 
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and the constraint in (13) reads Xlm = -^(^+) + ^{fi-) = -^i 

(n+) + = 1. (19) 

This pair of self-consistency equations, (18) and (19), is formally identical to that of 

the SU{2) AIM at UHF level [4]. In addition to the RHF solution described above 
{{n±) = i), a symmetry-broken solution is thus found to arise [4] whenever U > tt: 
physically the '+' and '— ' groups are split about w = 0, leading to a net lowering 
of the energy. The splitting of the levels is described most easily by introducing 
the 'moment' /i, writing (n±) = ^(1 T m) (which automatically satisfies (19)). The 
remaining self-consistency equation, (18), then reads 

/^=^tan-i(^^j, (20) 

from which it follows that symmetry-brokc^n solutions with /x 7^ arise in degenerate 
pairs; namely fi = +\i^o\ and fx = — |/io|, which we denote respectively as the 'A'-type 
and 'B'-type solutions. From (14) one sees that an 'A'-type solution has the energy 
levels at e± = ±5?7|/io|, while the 'S'-typc solution by contrast has e± = ^^U\iJ,o\. 

For a given type of solution {A or B) , the mean-field ground state is degenerate. 
The 2N mean-field impurity levels can be partitioned between two groups of N levels 
m Wn (= C^)) distinct ways. For example, one solution corresponds to allocating 
flavours m = 1 - N to the '— ' group and the remaining flavours m = {N + 1) - 2N 
to the group; a second solution arises when the two sets of flavour components 
are completely interchanged; and other solutions can obviously be obtained by 
interchanging only some of the flavour components between the two groups. Bearing 
this in mind, one can if desired enumerate all Wn possible degenerate MF states for 
any given A^, labelling them using the following notation. For each state, the 2N 
possible flavour components are written in a row, such that the first N belong to the 
' — ' group and the remaining N are those of the group. For example, in the case 
= 2 [5'C/(4)], the state (1 2 3 4) has ei = e_ = 62, and ej = e+ = £4. All other states 
can be generated by permutations of the fiavour components: e.g. the permutation 
(1 3 2 4) gives the state in which ei = e_ = €3, and €2 = e+ = €4. Note that some 
permutations lead to physically equivalent states e.g. (1 3 2 4) and (3 12 4) — since 
the ordering of levels within a group is irrelevant. To avoid this, one considers only 
permutations in which the two sets of N integers (corresponding to the '— ' and 
groups) are each in ascending numerical order; i.e. the state (1 3 2 4) is allowed, 
whereas (3 12 4) is not. In this way, every distinct MF state is counted once, and 
once only. 

The net MF Green fimction for flavour to, denoted here by G'^{uj), then follows 
by averaging the broken-symmetry, flavour-m Green functions of the Wn degenerate 
ground states: 

G^(a;) = :ji-^g;j,'^(c.) (21) 

where Q!^^^{uj) is the flavour-m Green function for MF solution a (= either A or B) 
and a particular permutation P of the flavour components between the e+ and e_ 
levels. The sum over P is taken over the set of permutations that enumerates all Wn 
physically distinct states, as described above. Equation (21) can then be simplified 
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by the following arguments. First, as seen explicitly in (15), the MF Gmi^) depends 
on Cm only. Therefore, on defining 

S±{uj) = — rTl^> (22) 

we obtain CJ^^{u}) = G±{ijj) if the permutation P allocates flavour component m to 
the '±' group, respectively. And by symmetry, when the Wn degenerate states are 
enumerated the flavour m. will appear Wn/2 times in each of the '+' and ' — ' groups 
{Wn is necessarily even for > 1). Combining these results, one obtains 

GlM = ^[e?M + ^?!!M] (23) 

for cither a = A or B (and obviously being independent of to, reflecting SU(4) 
rotational invariance). 

Hence, the net MF Green function is equivalently the average of the MF Green 
functions for the + and — groups, regardless of how the flavour components are 
actually allocated to those groups to form a particular MF solution. Moreover, 
since intercoverting A and B exchanges the energies of the two groups such that 
e±(w) = Q^iio), (23) becomes 

G^M = ^[^±M + elM]; (24) 

the statistically averaged MF Green function is thus equivalently the average of the 
'A'- and 'i?'-type Green functions for either given group (+ or — ), again regardless 
of how the flavour components are partitioned between the groups. In either case, on 
taking the imaginary part and using (22), the average UHF spectrum follows as 



A2 A2 
^0 _|_ ^0 



(25) 



.(u;-iC/|/zo|)2+A2 (c. + iC/|Mo|)2 + A2. 

consisting therefore of two superimposed Lorentzians ('Hubbard satellites') of half- 
width, half-maximum (HWHM) Aq. 

Note that when Ag is switched off, the MF spectrum (25) reduces correctly to 
that of the atomic limit, itself obtained from the imaginary part of (10). This is 
seen directly by taking the limit Aq ^ in (25): the self-consistency condition of 
(20) becomes n = sgn(/u) [and hence |/xo| = 1], while the Lorentzians become delta- 
functions at a; = i^C/l/xol = ^^U. And when U = 0, the non-interacting spectrum, 
(9), is trivially recovered by (25). 

The UHF approach thus captures both the atomic and non-interacting limits 
(sees. (3.1,3.2)), and in some sense bridges the two when U and Aq are both finite. But 
it does not of course capture the full physics of the model. Equation (25) is independent 
of N; in reality, the widths of the Hubbard satellites increase with N, since this opens 
up more relaxation channels which lead to enhanced lifetime broadening [10, 19] (see 
also section 4.1 below). Most importantly, the low-energy description provided by (25) 
is entirely wrong. The spectrum of the AIM is well known to possess an exponentially- 
narrow Kondo resonance straddling the Fermi level [1]. This many-body resonance 
is simply absent at the static MF level of description, and can only be captured by 
developing a more sophisticated approach. 

In the next section we describe how the local moment approach can be extended 
to the SU{2N) AIM. As explained previously in a number of works (see e.g. [19,44]), 
the technique overcomes the intrinsic limitations of static MF theory via inclusion. 
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within an inherent two-self-energy framework, of dynamical self-energy contributions 

which in physical terms embody tunneling between degenerate MF states; and which, 
in acting thereby to restore the symmetry broken at pure MF level, lead correctly to 
recovery of Fermi liquid behaviour on low-energy scales. 

To set up the LMA in practice, it is convenient to work with the 'capacitivcly- 
coupled' quantum dot picture of (4), replacing all m indices by And just as 

the MF Green function G°^{uj) ((23) or (24)) is independent of the MF solution with 
which one chooses to work, the same can readily be shown to hold within the LMA. 
We thus choose to work with the particular MF solution in which all a =t levels 
are allocated to the '— ' group, and all a =1 levels are allocated to the group. 
Physically, this solution corresponds to uniform charges on all sites i, but broken spin- 
symmetry; more precisely, when > the 'A'-type solution has an excess of t-spin 
electrons on every (equivalent) site, while the 'B'-type solution has an excess of J.- 
spins. Other partitionings of the flavour components {i, a) would of course correspond 
to different physical pictures - such as MF solutions that are spin-symmetric but of 
broken charge symmetry - but we reiterate that which is used in practice is irrelevant 
when calculating the full rotationally invariant Green function. Writing (23) and (24) 
in the chosen charge-symmetric, broken spin-symmetry picture, gives 

Glioj) =l[Gr^{io)+grM (26) 

= IK{^) + GFA^)] (27) 

(where a = either A or B), which will be employed in due course; with symmetries 

t/j^(cL') = Qf_^{uj), such that from (27) G^^{lo) is independent of spin cr (and as such 
spin-rotationally invariant [19,49]), as well as being independent of i (reflecting the 
equivalence of the N levels in the capacitively coupled quantum dot picture). 

5.^. Local moment approach 

The MF solutions above provide a natural starting point for perturbation theory within 
a two-self-energy description [19,44,49]. We take the unperturbed Hamiltonian to have 
the same symmetries as that of MF, i.e. to be of form 

^0 = E ^kcLcH. - E + ^ E (4../.. + h.c.) (28) 

for each oi a = A ov B, taking A = + and B = — for notational convenience when 
using the label a; and where x = U\iJ,\/2, with local moment (determined at 

post-MF level via the symmetry restoration condition intrinsic to the LMA [44], and 
specifled at the end of the section). The perturbation term corresponding to (28) is 
Hf = H-H§, viz. 

Hi = '^Uhi^hii + ^ E ^"-^i E^^ "^^)^i<^' (29) 

i i,j i,a 

and in precise parallel to the MF results of the previous section, the unperturbed 
Green functions for a cr-spin electron on orbital i are 

, , 1 
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Figure 1. (a) Separation of tlie LMA self-energy T,ic,{uj) into static plus 
dynamical terms. The arrowed lines represent the unperturbed Green function 
Qia{ijj), and the wavy line represents the Coulomb interaction U. (b) Dynamical 
self-energy approximation within the LMA. In both figures the prime on the 
summation indicates that the term with j = i and a' = a \s omitted. 



The corresponding perturbed Green functions can be obtained by diagrammatic 
perturbation theory in Hf [49]. One naturaUy obtains the Dyson equation 

Gt,{Lo) = i . = (31) 

which defines the two self-energies Sf^(a') (for a arising for a given a [our 

definition of Yif^iw) excluding the static first-order diagram (e + aax) coming from 
the third term in (29)]. The full, rotationally invariant impurity Green function then 
follows using directly analogous arguments to those of the previous section, such that 
(c/ (26)) 

(32) 



(i,cr)). The spin symmetry of the Hamiltonian (28) 

(33) 
(34) 



(again independently of m 
implies that 

while particle-hole symmetry implies 

G^cr{(^) = ~G^cri-^), 

thus simplifying the calculation of (32) in practice. 

Following [19,44,49], we now expand the S^s perturbatively in their respective 
s. For specificity we consider a = A only (and for clarity drop the A superscripts 
from here on unless indicated otherwise explicitly); the i?-type self-energy follows (if 
desired) from the symmetry (33). As in previous LMA papers {e.g. [19,44]), we choose 
to separate the self-energies into static (tj-independent) plus dynamical terms. The 
former is approximated by the sum over 'tadpole' diagrams in figure 1(a), which gives 



Si<T(w) 



^U{l + <j'\fi\) + ^M 



{N-^)U-^aU\fi\+J:M 



(35) 
(36) 
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Figure 2. Transformation of figure 1(b) using the fact that Qjai is independent 
of j. The shaded bubble represents the RPA spin polarization propagator, (40). 



with 



IMI 



■ tan 



(37) 



(and X — x/Aq). Generalising the approach of [19], we approximate the dynamical 
part of the self-energy for a flavour-m [= (i, a)] particle by the infinite series shown in 
figure 1(b). As we are working from the broken-symmetry states in which the Qja'S are 
independent of orbital index j, all N diagrams describing interactions with — cr-spin 
electrons are equivalent, as are all remaining — 1 diagrams describing interactions 
with cr-spin electrons. The sum over j, a' in figure 1(b) can thus be replaced by just 
two terms. Upon recasting the diagrams in terms of random-phase approximation 
(RPA) polarization propagators, we obtain the result shown in figure 2. 

Note that all explicit A^ dependence to the self-energies is contained in the 
coefficients of the two diagrams, the broken-symmetry Green functions from which 
the diagrams are constructed being independent of explicit factors of A^. In the SU{2) 
case (A^ = 1), only the first diagram on the right-hand side of figure 2(a) survives 
and the self-energy approximation reduces to that used in [19]. When A^ > 1, both 
diagrams contribute, although we find that the dominant contribution in the strongly 
correlated regime, t/ ^ 1, comes from the first diagram alone [see Appendix B]. 

Translating the diagrams in figure 2 using the Feynman rules gives (c/. [19]) 

S],.(u;) = Nl^l-'^ico) + {N- l)^mco) (38) 



where 



27ri 



The RPA polarization propagators H'^'^ (uj) appearing above are translated as 

(w) 



1 - [/ on" 

with the 'bare' polarization bubble 

1 
27ri 



dWl 0io-'(^l)^icr(^l -^)- 



(39) 



(40) 



(41) 
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For the flat-band AIM of interest, the latter can be obtained in closed form. The 

equations for '^H'^^'^ (uj) arc provided in the Appendix of [19], and those for ^'n"'"'{(j) 
are given in Appendix A of the present work. 

We note that the numerical calculation of the LMA dynamical self-energies is 
quite straightforward. In practice it is convenient to determine separately the real 
and imaginary parts of the various constituents listed above, the former obtained by 
Hilbert transformation of the latter. As the methods used are essentially identical to 
those of previous work [19], we do not discuss them further here. 

It remains finally to specify the condition for determination of the local moment 
(entering parametrically via x = U\ijl\/2 in the Hamiltonian (28)). This is achieved 
by 'symmetry-restoration', the central element of the LMA, corresponding physically 
to restoration of the symmetry which is broken at pure MF level. It amounts to a 
single self-consistency condition on the two self-energies precisely at the Fermi level, 
a; = 0; which, regardless of the specific diagrammatic approximation used for the self- 
energies, ensures that the rotationally invariant LMA Green function Gia{co) recovers 
Fermi- liquid behaviour on the lowest energy scales [44,49]. Since symmetry restoration 
is discussed extensively elsewhere (sec e.g. [44,45,49,52,58]) we do not elaborate 
further on it here; but simply note that in direct correspondence with earlier work, 
the symmetry restoration condition here is 

ti^{LJ = 0) = ti^{iO = 0). (42) 

In addition, analyticity of the RPA transverse spin polarization propagator requires 
[19] 

2x 

ImoI < ImI = < 1. (43) 

Equations (42) and (43) are readily solved numerically for any choice of the bare model 

parameters, and are suflicient to determine x uniquely in all cases. 

For later use, we also note that a combination of particle-hole symmetry [(34)] 
and t / i-spin symmetry [(33)] implies that (42) reduces to 

e + Ei<,(a; = 0)=0. (44) 

Hence, from (36) and (6), the symmetry restoration condition of (42) can be written 
alternatively as 

Si.(a; = 0) = |[/|A|. (45) 

4. Results 

We have analysed in detail the SU{2N) AIM within the LMA described above. 
Numerical evaluation of the impurity Green function is computationally very 
inexpensive, and we have examined its behaviour over a wide range of the bare model 
parameters. Here we describe our main results, focussing in particular on the strong- 
coupling limit, U = U/Aq ^ 1, where the model enters the Kondo regime in which 
electron correlations are key. 

4.1. Single-particle spectrum: overview 

Figure 3(a) shows representative single-particle spectra of the SU{4:) AIM in strong- 
coupling, calculated within the LMA as described in the previous section, for U = 20, 
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Figure 3. Representative single-particle spectra of the SU{2N) AIM: 5f7(4) 
results arc shown for U = 20, 40 and 50 (solid, dashed and dot-dashed lines, 
respectively), and the 5(7(2) spectrum for U = 20 is also shown as a dotted 
line, (a) shows the behaviour on high frequency scales lo ^ U , while (b) and (c) 
highlight the low-frequency Kondo resonance. 

40 and 50. The spectra are plotted as 7rAoD(w) vs w/Aq on a linear frequency scale, 
and at the Fermi level in particular are seen to satisfy correctly the dictates of the 
Freidel sum rule for all SU{2N) at half-filling, 7rAoL'(w = 0) = 1 [1,24]. The figure 
also shows for comparison the corresponding spectrum of the SU{2) AIM [19] for 
U = 20. 

Consider first the high-frequency behaviour seen in figure 3(a). The spectra all 
show clear Hubbard satellites arising on the scale Iwj ~ U/2, reflecting physically the 
single-particle excitations of the atomic limit (see section 3.2). The finite coupling 
to the host broadens the Hubbard satellites from the atomic-limit delta-functions in 
(10). This broadening is ©(Aq) in both cases, but the Hubbard satellites of the SU{4:) 
spectra are significantly broader than those of the SU{2) spectrum (moreover, both 
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are broader than the UHF spectrum in (25)). The origin of this effect, as alluded to 
at the end of section 3.3, is many-body scattering from the impurity site [19], which is 
not of course captured at pure MF level. The energy cost to excite particles between 
different flavour states on the impurity site {i.e. 'flip a spin' in the SU{2) case) is of 
order the Kondo scale. This is negligible on the energy scale of the Hubbard satellites, 
and hence the number of relaxation channels contributing to the lifetime broadening 
of the Hubbard satellites increases with A''. Indeed in the particle-hole symmetric limit 
considered here, it is straightforward to show from the LMA that the widths of the 
Hubbard satellites go as {N + l)Ao in the strong coupling limit, U ^ 1 [10,38]. 

Moving to lower energy scales, we see that the LMA indeed captures the SU{4:) 
Kondo physics, all spectra in figure 3(a) displaying a sharp Kondo resonance centred 
on the Fermi level u = 0. This behaviour is clearer on the much reduced linear 
frequency scale shown in flgure 3(c); but is seen most clearly by showing the spectra 
on a logarithmic frequency scale as shown in figure 3(b).§ The evident conclusion 
to be drawn from figure 3(b) is that the low-energy behaviour of the three SU{4:) 
spectra (solid lines) is universal [38,42]. Defining a characteristic Kondo scale wk by 
the HWHM of the Kondo resonance, it is quite clear that the three solid lines will 
scale onto each other when plotted on the reduced frequency scale lu/lok (provided 
naturally that u) itself does not approach non- universal scales of order |a;| ~ 0(Ao)); 
as considered in detail in section 4.3 below. Note further that the SU{2) spectrum 
shown as a dotted line in figure 3(b), itself known [19,45] to display universal scaling 
as a function of ui/lok, does not scale onto the universal SU{4:) behaviour: i.e. the 
SU{2) and S'C/(4) scaling spectra are distinct [42] (c/. [38]). To understand these 
results further, we now take a closer look at the low-energy physics of the model. 

4.2. Low-energy scale 

Spectral scaling reflects of course the single low-energy scale inherent to the problem 
in strong coupling {/ » 1, viz. the characteristic SU{2N) Kondo scale. This arises 
naturally within the LMA, in the imaginary part of the polarization propagator 
11"' (w) (see (40)), where it is manifest as a sharp resonance with maximum at 

1^ = l^m OC U)K [19]. 

That the LMA is exponentially small in [/ is a direct consequence of the well- 
known incipient divergence of the RPA ladder sum in (40) , together with self-consistent 
enforcement of symmetry restoration, (42). One can in fact extract the leading U 1 
behaviour of the scale analytically, following [19]. Details of the calculation are given 
in Appendix B, where we obtain the key result 



with c{N) a constant prefactor for a given N. The exponent agrees with numerics 
[38,40] and a slave-rotor mean-field theory [37], and it contains the same (exact) 
1/N dependence known from analytical studies of the SU{2N) model in the U ^ 00 
limit [1,10,27]. This 1/N factor in the exponent means of course that the Kondo 
scales for SU{2N) models with the same U ^ 1 but different N are vastly different, 
as indeed seen explicitly in figure 3. 

§ The small spectral 'glitch' (seen e.g. at ti; ~ 10~^ in the solid line of figure 3) is a well known 
artifact of the RPA ladder sum used within the LMA. While its effect on the spectrum is slight, it 
can be removed if so desired by the procedure described in [45]. 




(46) 
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Figure 4. The SU{2N) Kondo scale, plotted as ln{u!m/U) vs. U. The points are 
numerical results for the SU{2), SU{A) and 5C/(6) models (triangles, squares and 
circles, respectively), while the lines are fits to (46) for each group. 



To verify (46) itself, we have calculated numerically over a wide range of 
i7 > 1, for the SU{2), SU{A) and SU{<o) cases. The results are shown in figure 4, 
plotted as \n{uo^/U) vs. U . The solid lines are the predictions of (46), with the constant 
c{N) fit to the numerically determined points (and found in practice to depend rather 
slowly on N). The agreement between the numerics and (46) is excellent. 

^.3. Spectral scaling 

Having identified the Kondo scale Wm, we return to the issue of spectral scaling in more 
detail. Figure 5 shows explicitly the scaling collapse of the SU{4) Kondo resonance as a 
function of w/wm, for U — 30, 40 and 50 (all collapse perfectly to the scaling form over 
the uj/uia range shown). For comparison, the LMA SU{2) scaling spectrum [19,45] is 
also shown, dotted line. The figure highlights the clear A^-dependence of the SU{2N) 
scaling spectrum, alluded to in section 4.1. The SU{A) Kondo scaling resonance is 
seen to be somewhat narrower in form compared to its SU (2) counterpart (although 
it is of course broader on an 'absolute' w-scale, reflecting as discussed above that 
uJm{N = 4) > Wm(Af = 2) for given U). 

The scaling behaviour of the tails themselves is considered in figure 5(b). Here 
we show the universal forms of 7rAo-D(w) vs. uj/uJni for the SU{2), SU{4:) and SU{6) 
cases, plotted with logarithmic axes to highlight the A^-dependence of the spectral 
tails (the corresponding linear plot is shown in figure 5(c) for comparison). As N 
increases, the tails are seen to decay more rapidly, although the form of the decay is 
evidently similar in each case. 

This behaviour of the scaling spectrum, for |w|/wm ^ 1 and arbitrary A^, can in 
fact be obtained in closed form, following closely the approach of [45] . Details of the 
calculation are given in Appendix C: the final result is found to be 



7rAoi:'(w) 



1 1 4A^ + 1 



(47) 



2 [{^ln\Cj\)' + l (MlnM)%(4A^+l)2 
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Figure 5. Universal scaling spectra, (a) Scaling collapse of the SU{A) spectrum 
in strong coupling [0 2> 1), as a function of cj/ojm, for U = 30, 40 and 50 (solid, 
dashed, and dot-dashed lines respectively). Spectral collapse is essentially perfect 
on the energy scales shown. The dotted line shows the corresponding universal 
spectrum for the SU (2) model [19] . 

(b) Scaling spectra for SU{2), 5(7(4) and SU{6) (solid, dashed and dot-dashed 
lines) on a logarithmic scale, together with comparison (dotted lines) to the 
analytic forms in (47); 

(c) shows the decay of the spectral tails for 5(7(2), 5(7(4) and 5(7(6) on a linear 
scale. 
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Figure 6. Comparison between LMA and NRG scaling spectra. Solid and dashed 
lines show respectively the LMA 5(7(2) and SU{4:) scaling spectra, while the 
dotted lines are corresponding NRG results. The Kondo scale uj-k is defined by 
7rAoD(a;K) = 1/2. Note that the SU{4) spectra have been shifted upward by 0.4 
for clarity. 

where lj — lu/lo^. The asymptotic resuhs for iV = 1, 2 and 3 are shown m figure 5 as 
dotted fines, and are seen to agree excefiently with the numerical curves for \u)\ > 5. 

4.4- Comparison with NRG 

We conclude our discussion of the SU{2N) AIM by comparing the scaling spectra 
obtained above with essentially exact numerics from NRG [7,8,39,63,64], employing 
the full density matrix formulation of the method [63,64], and with the self-energy 
determined directly [65] . At the time of writing, available computational power renders 
it impossible to obtain accurate NRG results for models with N > 2. Nevertheless, 
the LMA scaling spectrum for = 1 is known already to agree very well with the 
NRG result [45], and here we show that the same is true for N — 2. As such, it is 
difficult to imagine the situation changing significantly for higher values of N. 

Figure 6 compares the LMA scaling spectra for the SU{2) and SU{A) models 
(solid and dashed lines respectively) with those of NRG (dotted) ; shown as a function 
of w/ujk with wk defined by ttAqD{lok) = 1/2- Note that both S'C/(4) spectra have 
been shifted upwards by 0.4 for clarity. 

The agreement between the LMA and NRG scaling spectra is striking. The tails 
of the spectra for |w|/wk ^ 1 are essentially indistinguishable, suggesting that the 
asymptotic form in (47) may be exact. While the LMA and NRG spectra begin to 
deviate when uj/ujk is reduced (in part due to the RPA 'glitch' in the LMA spectra), 
symmetry restoration ensures the correct spectral pinning required by the Friedel sum 
rule, viz. ttAqD{oj = 0) = 1, and hence the LMA and NRG spectra once again coincide 
as w ^ 0. 
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5. Away from SU{2N) symmetry 

The discussion of sections 2-4 has naturahy focussed on the fully-symmetric SU{2N) 
Hamiltonian of (4) . We now comment briefly on the effect of an anisotropic interaction 
between particles on the impurity. Taking (4) and replacing the coefflcient of the final 
term by U'/2, viz. 



U' v^' 

Y 2^ ninj,{A8) 



leads to a model of capacitively-couplcd SU (2) AIMs with distinct onsite and intersite 
interactions, U and U' respectively [40 42]. For U' ^ U, the Hamiltonian has SU{2)^ 
symmetry rather than the full SU{2N) symmetry of the parent model. We shall focus 
here on the case U' < U for simplicity. 

The consequence of explicit symmetry breaking in the Hamiltonian (48) is that 
for U' < U, the lowest energy broken symmetry mean- field solutions arc now only 
those with uniform charge on all sites i and broken spin-symmetry. On repeating the 
arguments of section 3.3 however, the structure of the MF Green function for U' < U 
remains precisely that given in (26,27). The LMA can then be set up in the same way 
as described in section 3.4, taking care now to distinguish between U and U' vertices 
in figure 1. Since each term in figure 1(b) involves repeated interactions between (i, a) 
and the same (j, a'), each term involves cither U or U' vertices, not a mixture of the 
two. It is then readily shown that the LMA self-energy of (38) generalises to 

^.Aoj) = U) + {N- 1) U') + U')] (49) 

where the U or U' vertices appearing in the various terms are labelled explicitly. 
For the reasons discussed in Appendix B, the contribution of the T,f^{ui; U') term in 
(49) is negligible when U ^ 1. The two remaining terms by contrast both involve the 
transverse spin-polarization propagator 11'^"'^ (w) , the resonance in which generates the 
low-energy Kondo scale within the LMA [see section 4.2]. Neither term can therefore 
be neglected a priori. And since each is calculated with different interaction vertices, 
there are now two distinct energy scales in the problem when U' ^ U . 

We have studied the ensuing U' < U behaviour in detail, both numerically and 
where possible analytically. The typical behaviour of the two energy scales is shown 
for the illustrative case of SU{4:) {N = 2) with [/ = 20 in figure 7(a). We denote by 

the scale that enters the first term in (49) (arising from spin-flip excitations on 
the same site), shown in figure 7(a), while the other scale (due to spin-flip excitations 
involving distinct sites) is labelled ui'^ and also shown in the figure. 

The two-scale behaviour is in good agreement with NRG results [40,41]. The low- 
energy scale (jjm deviates very little from its f/ = SU{2) behaviour [oc exp(— Trf/ /8), 
(46) with iV = 1] on initially increasing U' /U . Only when U — U' itself becomes of 
the order of the SU{A) Kondo scale does the scale rapidly cross over to the U' = U 
behaviour [cjm oc exp(— 7rl7/16)] discussed in section 4.2 [40]. The higher energy scale 

is typically of order U — U' , and arises in the NRG as a crossover scale from the 
SU{4:) local moment (LM) fixed point to the SU{2) LM fixed point [41]. And as 
U' U, figure 7(a) shows that lj'^ Wm, reflecting the natural existence of only a 
single scale at the 5i7(4) point U' = U. 

From our current perspective the most interesting physics is seen when the two 
scales are comparable; i.e. when U — U' 0{(jJm)- To express this more precisely, we 
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Figure 7. Results away from SU{2N) symmetry, (a) The behaviour of the two 
energy scales LJm/Ao (solid line) and w'^/Ao (dashed line), as a function of U' /U, 
for fixed U = 20 and N = 2. (b) The universal single particle spectrum ■kAqD{u)) 
calculated close to the SU (4)-symmetric point, with a = (U — U') /LOm.{U' = U) = 
5 (solid line). Also shown is the pure SU{2) scaling spectrum (dashed line), and 
the pure SU{A) scaling spectrum (dotted). 



define the quantity 

U -U' 



with uj„i{U' — U) the SU{A) Kondo scale arising when U' = U. For fixed N we 
find that when U ^ 1, the ratios lo'^/uj^ and Wm/'^m(t^' ~ U) ^ 7(a) are universal 
functions of a for any finite a. As a result, spectra for fixed a > show universal 
scaling collapse in the same way as seen previously for a = (section 4.3). 

We show a representative example of one of these scaling spectra in figure 7(b), 
for the case a — b (solid line). Note that in addition to the spectral 'glitch' at 
uj — tJm, there is now another at w = ui'^ ~ 25aJin- Both are a simple consequence of 
the RPA self-energy approximation (and can be removed if so desired [45]), although 
in this case they serve a useful purpose in marking out the low-energy scales of the 
model. Figure 7(b) shows in addition the universal SU{2) scaling spectrum as a 
function oiu/u)^ (dashed line), and the universal SU (4) spectrum also as a function of 
w/wni = w/[7(a)wni(f^' — U)] (dotted fine). Just as found in NRG studies [42], we see 
immediately that the U' < U scaling spectrum shows a clear crossover from universal 
SU{4) behaviour when \uj\ » w^, to universal SU{2) behaviour when |u;| <C w^, 
reflecting the different effective low-energy models that arise on the two distinct energy 
scales [40,42]. 

6. Conclusions 

We have considered here a local moment approach to single-particle dynamics of the 
orbitally degenerate SU{2N) Anderson model; focussing explicitly on the particle- 
hole symmetric case for arbitrary N, where the impurity is occupied by N electrons 
- a tangible example being the middle of the 2-electron Coulomb blockade valley in 
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a carbon nanotube quantum dot with S'C/(4) symmetry [21-24]. The resultant LMA 
is a natural generalisation of the technique developed previously for SU{2) quantum 
impurity models; and is seen to provide a rather good description of the dynamics 
of the SU{2N) Anderson model, recovering the correct exponential vanishing of the 
Kondo scale for J7 3> 1 and its associated universality, the broadening of the Hub- 
bard satellites by many-body scattering processes, and a very good description of the 
single-particle scaling spectrum as judged by direct comparison to NRG results for 
the SU{A) and SU{2) models. 



Appendix A. Closed form of °n<^'^(a;) 

The longitudinal spin polarization bubble '^11'^'^ (w) is calculated from (30) and (41). 
We find 



sgn(a;)7rAo Im (w) = — ^ v-; ■ ^v-; 2) 



where 



/((D) = tan (a) + i) — tan (a) — x) (A. 3) 

9{^) = In I 1 (A.4) 

/i((D) =(D(a)^+4). (A.5) 

Appendix B. Derivation of (46) 

The derivation of (46) for arbitrary N is closely analogous to that of the = 1 
case [19]. We begin by noting from (20) that the MF moment |/io| 1~ when 
f/ — > 00. Since stability of the RPA polarization propagators requires (43), symmetry 
restoration takes place at an x ~ U/2 and thus, from (37), \p,\ — »■ 1~ also. Hence the 
symmetry restoration condition, (45), becomes 

S.,(. = 0)^^>^|, (B.l) 



or, equivalently, 

NEffico = 0) + {N- l)S++(a; = 0) | (B.2) 
from (38). 

The first term on the left-hand side of (B.2) is simply N times the LMA self- 
energy for the SU (2) AIM, the asymptotic [/ » 1 behaviour of which has already 
been examined in detail [19]. The key is to note that the RPA n+~(a;) is dominated 
by a sharp resonance at u> = Wm, which as f/ — > 00 tends asymptotically to a delta- 
function, leading to [19] 

E+-(. = 0)^^-1^ In (-f). (B.3) 
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The second term on the left-hand side of (B.2) can be analysed in a similar fashion, 
starting from the closed form of "n"'""'"(ti>) in Appendix A, and using (40) to obtain the 
RPA longitudinal spin polarization propagator n++(a;). Crucially, in marked contrast 
to the n+~(w) discussed above, we find that no sharp resonance develops in the RPA 

n+"'"(cj) when U ^ 1. The numerically-determined 'Ef-^~^{uj = 0) thus turns out to be 

0{U^^), and hence in the limit U ^ 1, T,f^'^{uj = 0) can be neglected in comparison 
to the logarithmically-diverging T,f-^~{u> = 0) of (B.3). Combining (B.2) and (B.3) 
then leads to (46) straightforwardly. 

Appendix C. Derivation of (47) 

The tails of the SU {2N) scaling spectra can be obtained in closed form within the 
LMA. We define w = Lo/cOm and consider the behaviour of the LMA self-energies for 
fixed u) in the formal scaling limit 0. Let us rewrite (38) explicitly as 

^iA'^oj^) = NHI-^wuj^) + {N- 1)111^ {wuj^). (C.l) 

As described in Appendix B the imaginary part of the longitudinal spin polarization 
propagator, Imn'^'^(a;), contains only broad resonances at w ~ 0{U) rather than the 
emerging delta- function at w = in ImH"^'^ [lo). This means that Ef^(w) is non- 
universal at low frequencies and hence, in the limit — > for fixed a), we can write 
T.1^{u)UJ^) = Sf^''(0). But since Sf<^'(0) = 0{U-^) (as in Appendix B), it vanishes in 
the formal scahng limit U ^ oo. Hence (C.l) reduces simply to 

^,A^uj^) N^-'^icow^). (C.2) 

Moreover, {ujcu^) is precisely the function already analysed to obtain the 

scaling spectrum of the SU (2) AIM [45] . All the analysis of [45] can thus be generalised 
to the SU{2N) model quite straightforwardly. The upshot is that eqn. (5.1) of [45] 
must simply be multiplied by N (refiecting the additional factor of N in (C.2)) which, 
by analogy to eqn. (5.2) of [45], gives the full LMA scaling spectrum for the SU{2N) 
AIM as 

[Mln|^ + l|]Vl [ifln\c,-l\Y + [l+m0{u,-lf\ 

for w > (the corresponding negative-frequency spectrum follows from particle-hole 
symmetry, i.e. D{—w) = D{w)). In the limit [wj » 1, (47) results. 
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